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ABSTRACT 

In this article we extend the study performed in our previous article on the collapse 
of primordial objects. We here analyze the behavior of the physical parameters for 
clouds ranging from IO^Mq to IO^^Mq. We studied the dynamical evolution of these 
clouds in two ways: purely baryonic clouds and clouds with non-baryonic dark matter 
included. We start the calculations at the beginning of the recombination era, following 
the evolution of the structure until the collapse (that we defined as the time when the 
density contrast of the baryonic matter is greater than 10'’’). We analyze the behavior 
of the several physical parameters of the clouds (as, e.g., the density contrast and 
the velocities of the baryonic matter and the dark matter) as a function of time and 
radial position in the cloud. In this study all physical processes that are relevant to the 
dynamical evolution of the primordial clouds, as for example photon-drag (due to the 
cosmic background radiation), hydrogen molecular production, besides the expansion 
of the Universe, are included in the calculations. In particular we find that the clouds, 
with dark matter, collapse at higher redshift when we compare the results with the 
purely baryonic models. As a general result we find that the distribution of the non- 
baryonic dark matter is more concentrated than the baryonic one. It is important to 
stress that we do not take into account the putative virialization of the non-baryonic 
dark matter, we just follow the time and spatial evolution of the cloud solving its 
hydrodynamical equations. We studied also the role of the cooling-heating processes 
in the purely baryonic clouds. 
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1 INTRODUCTION 

The models of structure formation consider that all the 
structures of the Universe, such as galaxies and clusters of 
galaxies, were formed from the growth of small primordial 
density fluctuations present, initially, in an homogeneous 
and isotropic Universe, and that by gravitational instabili¬ 
ties resulted in the structures that we observe today. We are 
going to see below that we just start the study of formation 
and evolution of primordial clouds from density perturba¬ 
tions present at the beginning of the recombination era. 

The first galaxy formation models that included the ef¬ 
fects of gas dissipation were studied by Larson (e.g., 1969, 
1974, 1975, 1976). Recently, some authors used hydrody¬ 
namics codes including the gas dynamics (e.g., de Araujo 
1990, de Araujo & Opher 1988, 1989, 1991, Haiman et al 
1995, Thoul & Weinberg 1995, Tegmark et al 1997, Oliveira 
et al 1998, hereafter paper I, among others) in order to fol¬ 
low the behavior of the baryonic matter and to study the 
influence of several physical processes on the evolution of 
the primordial clouds. 

Gravity is generally believed to be the dominant force in 


the processes of structure formation and the dynamical stud¬ 
ies of these systems, based on observational data, reveal that 
the constituents of the Universe include not only luminous 
matter, but also a far greater amount of dark matter which 
influence the luminous matter through its gravitational field. 
For systems such as galaxies and clusters, hydrodynamical 
effects are important in the process of their formation and 
thus, hydrodynamics codes which include dark matter are 
then the most accurate manner to follow the cosmological 
evolution of the primordial clouds (see, e.g., Cen et al 1990, 
Cen 1992, Cen et al 1993). 

On the other hand, there are several non-gravitational 
processes that are important to the evolution of primordial 
clouds at high redshift, as the photon-drag due to the cosmic 
background radiation, the recombination processes and the 
hydrogen molecular formation. These physical processes are 
important during and after the recombination era. In gen¬ 
eral, the models that use hydrodynamical numerical codes, 
found elsewhere, do not initiate the calculations at the be¬ 
ginning of the recombination era, when the perturbations 
of the baryonic matter can effectively begin to be signifi¬ 
cantly amplified, whether they are Jeans unstable, and do 
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not include the above mentioned physical processes. 

In the present article we used the set of hydrodynamic 
equations of the paper I to analyze the behavior of the phys¬ 
ical parameters using first only purely baryonic clouds. We 
then compare the evolution of the physical parameters in¬ 
cluding non-baryonic dark matter. We start all the calcu¬ 
lations present here at the beginning of the recombination 
era and take into account all the relevant physical processes 
present during and after the recombination era, besides the 
expansion of the Universe. 

In §2 we describe the basic equations, in §3 we discuss 
our results and finally in §4 we present our conclusions. 


Compton cooling (see, e.g., Peebles 1968), is the cool¬ 
ing by molecular hydrogen (see, e.g., Lepp & Shull 1983) 
and La is the Lyman-a cooling (see, e.g., Calberg 1981). 
The temperatures in our calculations are never greater than 
~ 10, OOOA' and thus, the cooling-heating mechanisms con¬ 
sidered here are the dominant mechanisms for the thermal 
history of the baryonic matter contained inside the clouds. 

In order to consider the DM (non-baryonic dark mat¬ 
ter) models, we included the following equations: 


dt 
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2 EQUATIONS 


We consider spherically symmetric density perturbations 
which produces clouds of baryonic matter and dark mat¬ 
ter with densities greater than the density of the Universe. 
We consider the baryonic matter and dark matter as two 
fluids coupled by gravity. 

The equations for a purely baryonic matter are de¬ 
scribed in greater detail in paper I, and so, here we only 
summarize them, namely: 
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the various terms appearing in the Eqs.(5)-(8) are analo¬ 
gous to those found in the hydrodynamical equations to the 
baryonic matter. 

We consider that the primordial clouds evolve within a 
medium of density (the density of the Universe) and so 
we take into account this fact in our study; see paper I to 
the hydrodynamical equations of the Universe. 

The initial perturbation is taken to be in the form of a 
power law spectrum as considered, for example, by Gott & 
Rees (1975) 



(B, -B,) 

Pe ’-H (1 - x) 


^1+7 

mp J 


( 4 ) 



(1 -|- Zrec) 


( 9 ) 


As in paper I, the Eqs.(l) —(4) are, respectively, the 
equation of conservation of mass, conservation of momen¬ 
tum, energy and ionization balance. Where p is the mass 
density of the cloud, r is the radial coordinate, v is the ve¬ 
locity. The velocity is written as u = Vn + v^j, where 
is the peculiar velocity of the cloud and = Hr (where 
77 = R{t)/R{t) is the Hubble parameter, 77(7) is the scale 
factor and 77(7) its time derivative). The pressure of the 
cloud is P = kNpTmil + x) (where k is the Boltzmann 
constant, N is the Avogadro’s number, Tm is the tempera¬ 
ture of the matter, and x is the degree of ionization), is 
the Thomson cross section, b = do-g^/c (where Ugg is the 
Stefan-Boltzmann constant), Tr is the temperature of the 
radiation, rrip is the mass of the proton, c is the velocity of 
light. 

The gravitational potential is Vp = GM(r)/r^, with 
G the gravitational constant. The term A7(r) is the total 
mass (i.e. the baryonic and non-baryonic dark matter) con¬ 
tained inside the radius r of the cloud. We here take the 
total amount of baryonic matter and dark matter constant 
during the evolution of the clouds. 

The cooling-heating processes are included in the cool¬ 
ing function L. This term is the summation of four mecha¬ 
nisms; L = +La- Where, Lj^ is the cooling due 

to recombination (see, e.g., Schwarz et al 1972), is the 


with M the mass of the cloud, Mo the reference mass, n 
the spectral index and Zrec the redshift at the beginning of 
the recombination era. It is important to stress that M is 
the total mass of the clouds, i.e., M — Mg + M^^ (where 
“B” refers to the baryonic matter and the “DM” refers to 
the non-baryonic dark matter). The normalization of the 
spectrum of perturbations given above was made to agree 
with the COBE data (see e.g. de Araujo & Opher 1994). 

We consider that initially the density profile of the cloud 
is “top hat” like, and we begin the calculations when the ra¬ 
diation temperature is ~ 4, 00077 {z ~ 1500) where the ion¬ 
ization degree begins to be significantly lower than one. We 
use n_B = 0.1 (the baryonic density parameter) and different 
values for Qdm (the non-baryonic dark matter density pa¬ 
rameter), as well as different values for the Hubble constant, 
h (the Hubble constant in units of lOOfcm s~^Mpc~^). We 
consider that the clouds initially expand with a field velocity 
given by the Hubble flow. 

We follow the evolution of the cloud until its collapse, 
that we consider as the time, after the recombination era, 
when the density contrast of the cloud is greater than 10^. 
We used in our calculations 1,200 spherical shells for the 
baryonic matter and 1,200 shells for the dark matter. Mod¬ 
els with 6,000 shells for both components showed no consid¬ 
erable difference in the results. The complete set of partial 
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differential equations has been written in the form of differ¬ 
ence scheme by one of us (Miranda 1998) according to the 
basic method described by Richtmyer and Morton (1967). 

3 CALCULATIONS AND DISCUSSION 

The scale of masses studied ranges from lO^M© to 10'^®M©. 
We show, in particular, three points in the cloud; the first 
shell, the middle shell of the scheme, and the last shell as a 
function of time. In order to analyze the radial behavior, we 
show five different times for the output of the data: L - time 
at the recombination era, t 2 = (ts — ti)/2, ts = (tf — ti)/2, 
t 4 = {tf — tp)/2 and tf, where tf is the time of the collapse 
of the core. 

Also, as already mentioned in the paper I, due to the 
expansion of the Universe, even the clouds that are already 
Jeans unstable and eventually collapse undergo, initially, an 
expansion phase before collapsing (see paper I for more de¬ 
tails). 

We divided our paper in two parts. Firstly we stud¬ 
ied purely baryonic clouds and we analyzed the behavior of 
the variables along the evolution of time, and in different 
points inside the cloud. We also studied the influence of the 
cooling-heating processes. To do this we studied individu¬ 
ally the effects of all the four cooling-heating processes and 
the photon drag (Tables 2—4). The results are described in 
subsection 3.1. 

Secondly, we included DM and analyzed the behavior 
of the same variables in the same way, except the analysis of 
the cooling-heating processes. The results of baryonic plus 
DM model are described in subsection 3.2. 

3.1 Purely baryonic clouds 

We now show the results of the purely baryonic models. We 
ran our models for masses ranging from lO’^M© to 10^®M©. 
Here we ran the models with Mo = 4 x lO^^M©, h = 1, 
= 0.1 and n = — 1. The normalization of the spectrum 
used here is consistent with ag (see, e.g., de Araujo & Opher 
1994 for a discussion). We show these results in the Figs.1-6. 
It is worth mentioning that a spectrum such as J oc M~^t^ 
is consistent with observations for 1 < r/{h~^Mpc) < 100 
(see for example Kashlinsky & Jones 1991). 

Let us first have a look at the evolution of the density 
contrast, in the Fig.l, and also at the evolution of the den¬ 
sity, in the Fig.2. We see that masses within the interval 
lO^M© to 10^°M© have a partial collapse, i.e. our collapse 
criterion is satisfied only to the internal shells. In particular, 
for the masses lO^M© to 10®M© the density contrast has a 
slightly oscillation in the external shells and, these shells do 
not reach a density contrast of 10“^. 

We also analyzed the morphological behavior of the 
density profile (Fig. 3), in order to see how it is modified 
as a function of time. For the masses ranging from lO^M© 
to 10^®M©, the collapse starts in the internal shells, fol¬ 
lowed by the outermost shells, as discussed above. For the 
masses greater than 10^® M©, there is no collapse (at least for 
the normalization of the spectrum of perturbations studied 
here). Along the radius the initial “top hat” profile changes 
rapidly with time for masses smaller than 10^®M©. This is 
related to the fact that masses lower than 10^®M© have a 
pressure gradient that modify the initial density profile. This 


shows that to assume that the density profile is “top hat” 
throughout the cloud evolution is not a good assumption, as 
considered, e.g., by de Araujo & Opher (1989, 1991, 1994). 
We are going to see later on that under certain circumstances 
the use of a “top-hat” density profile throughout the calcu¬ 
lations, as assumed by de Araujo & Opher, can provide us 
interesting informations. 

Let us now have a look at the temperature, pressure 
and molecular density formation (Fig.4). These variables fall 
with time for all masses because of the cooling and the ini¬ 
tial expansion of the clouds. But when the cloud collapses, 
the temperature increases greatly as expected. An analy¬ 
sis of the molecular density versus temperature is shown in 
Fig.5. It shows that the molecular density increases when 
the temperature is within the interval lOOA < T < SOOA. 

As already explained in paper I, the heat conduction 
and the Bremsstrahlung are not taken into account. It is 
worth stressing that in our study we do not take, as initial 
condition, a virialized object as considered by some authors. 
Our calculations start at the recombination era when the 
temperature of the baryonic matter inside the perturbation 
has the same value as that of the radiation. For clouds that 
collapse at high redshifts, in general, the principal mecha¬ 
nism that acts in the evolution and collapse of the perturba¬ 
tions is the Compton heating-cooling (see results presented 
in Tables 5-8 in paper I). Thus, during the expansion phase 
of the perturbations, when the Compton heating-cooling is 
efficient, the temperature of the baryonic matter inside the 
perturbation is almost the same temperature of the radia¬ 
tion. 

At the redshift collapse of the clouds, the temperature 
of the matter increases by a factor 4 or 5 in relation to its 
correspondent value at the turn-around. Thus, depending on 
the power spectrum and on the mass of the cloud, the shocks 
and the collapse itself are not able to rise the temperature 
of the cloud to ICK (see, e.g.. Table 5 and the discussion 
in the Paper I). 

Certainly, these results are dependent on the particu¬ 
lar way in which we define the collapse of our models, that 
is, we consider that an object collapsed when the density 
contrast is greater than 10^. The results also depend on the 
thermal history of the clouds (that depends on the physi¬ 
cal processes here considered). As a result no strong shocks 
occur that could eventually heat the clouds to the virial 
temperature. If we follow the calculations to higher values 
of density than the established by our collapse criterion, 
the temperature will rise above the values here obtained. It 
is possible that to density contrasts higher than 10® strong 
shocks occur and they would probably heat the gas to values 
of temperature > lO^K, as a result other physical processes, 
as e.g. Bremsstrahlung, would be important to the evolution 
of the clouds. 

In Fig.6 we show the evolution, with the time, of the 
radii of the shells chosen to show the evolution of the clouds. 
It is possible to see the instant when the turn-around occurs 
for these shells. 

In Table 1 we present a comparison between our re¬ 
sults and the results by de Araujo (1990) and de Araujo & 
Opher (1994) for purely baryonic clouds. In this table we 
also present calculations using the normalization consistent 
with (Ts = 1 (see, e.g., de Araujo & Opher 1994), thereby 
we have to use Mo = 4 x 10^^ M© in Eq.(9). In fact this 
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normalization is more realistic since structures as large as 
M — may be formed. 

We note that all masses have collapse redshifts similar 
to those obtained by de Araujo and de Araujo & Opher using 
both normalizations, this is a very interesting result. If one 
would like only to have an idea of the collapse epoch, the use 
of a “top hat” profile throughout all the cloud evolution is 
a good approximation (see, e.g., de Araujo & Opher 1994). 

As mentioned in paper I the clouds stop expanding 
when (dp/p) > 4.6, therefore above the linear regime. As 
seen in paper I, we also obtain in the present calculations 
that the clouds stop expanding when (5p/p) > 4.6. This 
very fact explains why, even in the non-linear regime, some 
clouds studied here neither collapses nor stops expanding. 

It is worth stressing that the collapse redshift calculated 
here gives us only an idea of the epoch in which the collapse 
occurs, due to the fact that the formation of a galaxy, for 
example, is a process much more complex. In fact, what we 
do, it is a study of formation of a proto-object. The present 
study is of interest due to the fact that we show how a 
density perturbation could give rise to a proto-object and 
also we show which physical processes are important during 
and after the recombination era. 

We now study the influence of the cooling-heating pro¬ 
cesses on the evolution of clouds taking, as representative ex¬ 
amples, clouds with the masses lOM q, 10®M© and 
(Tables 2—4). To do this, we disregard, one at a time, all of 
the four cooling-heating processes and also the photon drag. 
The results are consistent with paper I, i.e., these very phys¬ 
ical processes are important on the evolution of the clouds 
determining, for example, the collapse epoch and the ther¬ 
mal history of the clouds. 

The analysis of the cooling mechanisms show that when 
we disregard the Compton cooling-heating, the collapses oc¬ 
cur earlier. This is due to the fact that the Compton mech¬ 
anism heats the cloud when the temperature of the matter 
is lower than the temperature of the radiation, in this way 
during the recombination era the Compton cooling-heating 
maintains the temperature of the cloud near to the temper¬ 
ature of the cosmic background radiation. Due to the expan¬ 
sion of the Universe the clouds firstly expand and whether 
the Compton heating was not present the clouds would ex¬ 
pand almost adiabatically, as a result their temperatures 
would get smaller, and so their pressures, thereby the clouds 
would stop expanding earlier. 

When we disregard the photon drag the clouds also col¬ 
lapse earlier. As we show in paper I the photon drag acts dis¬ 
tinctly on the smaller masses (< 10® M©) and on the greater 
masses. Here it acts against the collapse, as a brake, but its 
action is not so strong like it was for the smaller masses. 

The cooling by molecular hydrogen does not influence 
the collapse redshift, its influence is on the thermal history 
of the cloud. The other physical processes here studied have 
a behavior similar to that described in paper I. 

We studied also the behavior of the collapse of the 
clouds when we change h and fl. The results are in Ta¬ 
bles 5 and 6. As a general conclusion it can be mentioned 
that the collapse redshift, the maximum mass of the objects 
formed, and also the thermal history of the clouds depend 
strongly on h and fl. When we change, for example, from 
0.1 to 0.2 we have collapses occurring earlier. On the other 
hand, when we change h from 1 to 0.5 we have the collapses 


Table 1. The redshifts of collapse 



X 

II 

lOil’M© 

Mo 

= 1015M( 

M/Mq 

z* 

2^* 

Zc. 

2^ 

10" 

437 

- 

102 

108 

10® 

238 

- 

49 

49 

10® 

153 

- 

20 

21 

lOlo 

77 

75 

6.5 

7 

loll 

34 

34 

NC 

- 

1012 

12 

13 

NC 

- 

lOi® 

2.1 

3.6 

NC 

- 


We compare here the redshifts of collapse of the 
present article (z*), with those by de Araujo (1990) 

(zjf) and de Araujo & Opher (1994) (zjf*). We 
use n = —1, h = 1.0 and Qg = 0.1. We stud¬ 
ied two normalizations, namely, Mo = 10^® Mq and 
Mo = 4 X 10^’^Mq. NC stands for “no collapse”. 

occurring latter, as expected. 

3.2 Baryonic plus DM clouds 

We now discuss the results of our models including the pres¬ 
ence of the dark matter. We ran our models for masses rang¬ 
ing from IO^Mq to 10^® M©. Here we ran initially all models 
with Mo = 4 X lO^^M©, h = 0.5, = 0.1, = 0.5 and 

n = — 1, and we show these results in the Figs.7—13. In order 
to compare we also ran models with different combinations 
of h and , as well as, a different normalization of the 
spectrum of density perturbations. 

The behavior of the density is shown in Fig.7 (time evo¬ 
lution) and Fig.8 (spatial evolution). For the normalization 
of the spectrum of perturbations used here, the maximum 
mass that collapses is 10^'* M©. 

It is worth mentioning that the explanation of why there 
is no collapse for some clouds, as could be expected, is due 
to the fact, as already mentioned, that it is necessary to 
have (dp/p) > 4.6 to the cloud stops expanding and starts 
its collapse. 

We see, from Fig.7, rigth side, that all masses have 
partial collapse for the baryonic component, the outermost 
shells expand with the Universe. For the DM component 
(left side of the Fig.7), this partial collapse occurs only for 
the cloud of 10^^ M©. 

In Fig.8 we see that the top hat profile is destroyed as 
time goes on, particularly to the baryonic matter. In the 
DM component (left side of Fig.8) the top hat profile is 
maintained almost all the time during its evolution, but it 
is destroyed rapidly when the collapse is more evolved, in 
particular for the cases in which we have M > lO'^^M©. 

The DM profile is maintained due to the fact that it 
is pressure free, but when the gravitational influence of the 
baryonic matter, that is no pressure free, begins to be impor¬ 
tant the top hat profile of the DM component is destroyed. 

The behavior of the velocity profile is shown in Fig.9 
(time evolution) and Fig.10 (spatial evolution). In the bary¬ 
onic matter component (right side of Fig.lO), the profile of 
the velocity begins linear but it changes rapidly with time, 
just like in the purely baryonic model. The velocity profile of 
the DM component (left side of Fig. 12) is maintained linear 
almost all the time of the cloud evolution, only for clouds 
of M > 10^® M© there is a slightly modification in such a 
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Table 2. Values of the variables for M = IO^Mq for purely baryonic clouds. We used h = 1.0, 
Qb = 0 . 1 , Mo = 1 O^®M 0 and n = — 1 . 


Variable 

with all 

processes 

without 

La 

without 

La 

without 

without 

Let 

without 
photon drag 

2c 

103 

103 

105 

103 

103 

132 

(cm-3) 

10 ® 

10 ® 

1 

- 

10 ® 

300 

T{K) 

4 X 10® 

4 X 10® 

9 X 10^ 

8 X 10^ 

4 X 10® 

5 X 10® 

p{gcm 

3 X 10-1® 

3 X 10-1® 

3 X 10-1® 

3 X 10-1® 

3 X 10-1® 

5 X 10-1® 


Table 3. Values of the variables for M = 10® Mq for purely baryonic clouds. We used h = 1.0, 
Qb = 0 . 1 , Mo = 1 O^®M 0 and n = — 1 . 


Variable 

with all 

processes 

without 

La 

without 

La 

without 

^2 

without 

La 

without 
photon drag 

2c 

20 

20 

228 

20 

20 

28 


3 

3 

10-1 

- 

3 

1 

T{K) 

5 X 10® 

5 X 10® 

4 X 10® 

1 X 10"“ 

6 X 10® 

5 X 10® 

p{gcm 3) 

3 X 10-®i 

3 X 10-®i 

3 X 10-®i 

3 X 10-®i 

3 X 10-®i 

5 X 10-®i 


Table 4. Values of the variables for M = 10^^ M 0 for purely baryonic clouds. We used 


h = 1.0, LIb 

Variable 

= 0.1, Mo = 

with all 

processes 

10^^ Mq and 

without 

n = —1. 

without 

without 

without 

La 

without 
photon drag 

2c 

NC 

NC 

0.14 

NC 

NC 

2.6 


1 X 10-® 

3 X 10-® 

1 X 10-® 

- 

3 X 10-® 

3 X 10-® 

T{K) 

10 

50 

10 

10 

50 

4 X 10® 

p{gcm 

3 X 10-®® 

1 X 10-®® 

1 X 10-®® 

1 X 10-®® 

1 X 10-®® 

3 X 10-®"“ 


Table 5. Values of the variables for h = 1.0, Qb = 0.2, 
Mo = VC^Mq and n = — 1 for purely baryonic clouds. 


M/Mq 

Zc 

ng^{cm 3) 

T(K) 

p{gcm 3) 

iC 

110 

20 

4 X 10® 

3 X 10-®° 

10 ® 

24 

1 

6 X 10® 

1 X 10-®° 

loll 

2 

2 X 10-® 

3 X 10® 

3 X 10-®-! 

Table 6. 

Values of the variables 

for h = 0.5, 

Ub = 0.1, 

Mo = lOi® 

Mq and 

n = — 1 for purely baryonic clouds. 

M/Mq 

Zc 

"•Ha (cm 

T{K) 

p(gcm 3) 

10^ 

93 

300 

5 X 10® 

5 X 10-1® 

10® 

18 

500 

6 X 10® 

3 X 10-®® 

loll 

behavior. 

NC 

10-10 

10 

3 X 10-®° 


The temperature, pressure and molecular density for¬ 
mation of the baryonic matter component is shown in Fig.11. 
As can be seen in Figs.11-12, the production of the H 2 
molecules is efficient for all masses, once reached the tem¬ 
perature, of the baryonic matter, in the range lOOK < T < 
SOOTf. Near Zo the H 2 production increases, but for the ex¬ 
ternal shells (long dashed line) it is always low, due to the 
fact that the density of the cloud is also low in such shells. 

In Fig. 13 we show the behavior of the radii of the shells 
for the DM component and for the baryonic matter com¬ 
ponent (left side and right side of the figure, respectively). 


Comparing the two matter distribution it can be seen that 
the DM matter component is more concentrated than the 
baryonic matter component. 

At this point it is worth stressing that in our calcu¬ 
lations we perform a purely hydrodynamical study of the 
evolution of the clouds, we do not consider that neither the 
DM nor the baryonic matter undergo, for example, virial- 
ization. 

Due to the fact that there is no virialization for the 
DM, and that it is pressure free, it does not either undergo 
the photon-drag effect, and thus the DM gets more concen¬ 
trated. 

Its well known, and widely accepted, that the spiral 
galaxies have a halo with at least part of it made of non- 
baryonic dark matter (see, e.g., Persic & Salucci 1990). In 
this case the MbO/MbmO decreases with increasing r, 
contrary to the results present here, although we here are 
very far to have a complete model for a galaxy. Virializa¬ 
tion and star formation, for example, will certainly modify 
strongly the history of the proto-objets here studied. The 
standard scenario of galaxy formation outlined by, for ex¬ 
ample, White & Rees (1978) explains what could happen 
with the relative distribution of dark and baryonic matter. 
The dark matter forms a halo that collapses, and the gas 
shock heats to the virial temperature and is pressure sup¬ 
ported. The gas pressure then decreases through radiative 
cooling and then sinks leaving behind a halo of dark matter. 

In elliptical galaxies the distribution of dark matter is 
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not well determined as in the spiral galaxies. Based on the 
X-ray observations of hot gas of ellipticals, probably the best 
source of information in this case, some authors argued that 
Mbjjm ^ (where BDM is the baryonic dark matter), just 
like in the spiral galaxies. It is worth stressing, however, that 
these studies assume that the hot gas is isothermal, but usu¬ 
ally the temperature profile is poorly known (see, e.g., Carr 
1994). Whether the non-baryonic dark matter and the bary¬ 
onic dark matter have the same distribution, and whether 
the ellipticals follow the same DM distribution followed by 
the spirals, again there is no accordance between our results 
and the distribution of DM inferred for the ellipticals. On 
the other hand, we can argue that we are very far to have 
a complete model for an elliptical galaxy, to justify the lack 
of accordance. 

From X-ray observations it was shown that the cluster 
of galaxies A665, for example, presents a high concentration 
of (baryonic) dark matter at its center (see, e.g., Hughes 
& Tanaka 1992). Similar high concentration of dark matter 
at the center of clusters was derived for the Perseus (Eyles 
1992) and Coma (Briel, Henry & Bohringer 1992) among 
other clusters. Some authors argued that this dark matter 
should be dissipative and therefore baryonic, but it could 
well be, at least part of it, of non-baryonic nature. In this 
case the non-baryonic dark matter to baryonic matter distri¬ 
bution found in our study for structures with M > 10^^ M© 
would be in accordance with the distribution of DM inferred 
from the X-ray observations of clusters of galaxies, whether 
the non-baryonic dark matter follows a distribution similar 
to that followed by the baryonic dark matter. 

We also ran models with other combinations of the pa¬ 
rameters n, h, Udm and Mo to see how, for example, Zc and 
the maximum mass that collapses are modified. We show 
these results in Tables 7—9. 

In Table 7 we show the collapse redshift, Zc, for h = 1 
and 0.5 for the normalization Mo = 4 x IO^^Mq, as well 
as for the normalization Mo = 10'^®M©. For these cases we 
assume a Universe with Q.b =0.1 and Q.dm = 0.5. 

For the normalization Mo = 4 x lO^^M© we obtain no 
significant differences for Zc when we change h = 1 to 0.5. 
Also, the maximum mass for both cases is Mmax = 10^"^M©. 

For the normalization Mo = lO'^^M© we also obtain no 
significant differences for Zc when we change h = 1 to 0.5, 
for all clouds. 

Comparing now the results for the two normalizations, 
we obtain that the values of Zc are higher for Mo = 4 x 
lO^^M©. These results can be understood due to the fact 
that the normalization Mo = 4 x 10^^ M© have a greater ini¬ 
tial density contrast when we compare with Mo = 10^®M© 
and thus, for the same values of h and fl, the normalization 
Mo = 4 X 10^^ M© produces collapse of the clouds before 
those with the normalization Mo = 10^®M©. 

In Table 8 we compare our purely baryonic models and 
our baryonic plus DM models with the models by de Araujo 
& Opher (1991, 1994). For the purely baryonic cases, and 
in particular for the range of masses studied, there is no 
significant difference between the results for Zc present here 
as compared with de Araujo & Opher’s results. For the cases 
with non-baryonic dark matter we obtain values for Zc larger 
than that obtained by de Araujo & Opher. Comparing now 
only the results of our studies, we obtain that the inclusion 
of dark matter makes the collapses to occur earlier. 


Table 7. Models with baryonic plus dark matter 


Mo = 4 X IQi’^M© Mo = 10i®M© 


M/M© 

Zc 

Zc 


Zc 

Zc 

(h 

= 0.5) 

{h=l) 

(h 

= 0.5) 

(h=l) 

107 

515 

518 


131 

132 

10® 

332 

335 


66 

67 

10® 

191 

193 


31 

32 

lOlO 

100 

101 


14 

14 

loll 

49 

50 


5.7 

5.7 

lOi® 

22 

23 


1.6 

1.6 

1013 

9.4 

9.5 


NC 

NC 

10l3 

2.5 

2.2 


NC 

NC 

lOi® 

NC 

NC 


NC 

NC 

Collapse redshifts of baryonic 

plus 

DM models with 

Q.B = 0.1, U 

DM ~ ' 

0.5 and n = 

= -1. 




Table 8. 

Collapse redshifts of the models 



only 

baryonic 

baryonic plus DM 

M/M© 


Zc 

94 


•2c 

91 

10^ 

102 

108 

132 

113 

10® 

49 

49 

67 

56 

10® 

20 

21 

32 

- 

lOlO 

6.5 

7 

14 

- 

loll 

NC 

NC 

NC 

NC 


Comparison of the collapse redshifts of only bary¬ 
onic models (zc^), baryonic plus DM models 
(zcbm ) Araujo & Opher’s (1991, 1994) (zc^^ 

and Zc ) models with Mo = 10^®M©, n = —1, 

h = 1, = 0.1, = 0-5. 


Table 9. Comparison of the collapse redshifts of baryonic 
plus DM models with Mo = 4 X 10^’^M©, h = 0.5 and 

Qg = 0 . 1 . 


M/Mq 

n = —1 

{Qbm = 0-5) 

n = —1 
i^DM = 0-8) 

n = +1 
i^DM = 0 

10^ 

515 

520 

1471 

10® 

332 

336 

1435 

10® 

191 

194 

1366 

lOi® 

100 

102 

1235 

loll 

49 

50 

997 

lOi® 

22 

23 

635 

10l3 

9.4 

9.9 

270 

10i4 

2.5 

2.8 

77 

lOi® 

NC 

NC 

14 


From Table 9 we compare the results of the models when 
we change the slope of the spectrum (n) and the DM content 
Comparing the values of Zc when we use ~ 0 . 8 , 
instead of ~ 0 . 5 , we obtain no significant differences 
in the results. Also, the maximum mass in both cases is 
Mmax = 10 ^"' M©. Comparing now Zc, for ~ 0 . 5 , with 
n = — 1 and n = 4-1 we obtain very different values. The Zc 
values are very high for all masses studied when we adopt 
n = 4-1; in fact a power spectrum with such a spectral index 
produces excessive power. We have for this case collapses of 
masses of up to Mmax = 10^®M©. 
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4 CONCLUSIONS 

In the present investigation, we studied the evolution of 
density perturbations with the masses ranging from 10^ — 
IC^Mq, taking into account a series of physical processes 
which are present during and after the recombination era. In 
order to perform such a study, a spherical Lagrangian hydro- 
dynamical code was written, which made possible to follow 
the spatial and time evolution of the density perturbations 
with purely baryonic and baryonic plus non baryonic dark 
matter clouds. 

Our main conclusions for purely baryonic models are 
the following: 

a) Clouds with masses that range from —10^° M© have 
a partial collapse, for Mo ~ 

b) Clouds having masses > Mq do not have enough 
time to collapse, for Mo = 10^®M©; 

c) If one changes the perturbation spectrum we have col¬ 
lapse for masses up to lO^^M©; 

d) The Compton cooling-heating plays a key role in the 
collapse of the clouds; 

e) the photon-drag delays the collapse and 

f) The H 2 molecules are, in general, very important to the 
thermal history of the cloud. 

Our main conclusion for the baryonic plus DM models 
are the following: 

a) We have earlier collapse redshifts always, as compared 
with the purely baryonic models, independent of the combi¬ 
nation of the values of n, h, Udm and Mo that we used; 

b) In the collapse of the DM the cooling-heating mech¬ 
anisms and the photon drag, which are very important for 
the collapse of the baryonic matter, cannot directly alter the 
collapse of this component and thus the non-baryonic mat¬ 
ter is rapidly concentrated in the innermost region of the 
clouds; 

c) The observations suggest that in clusters of galaxies the 
dark matter is more concentrated than the baryonic matter. 
This result is in agreement with our results if the dark mat¬ 
ter is of non-baryonic origin. 
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FIGURE CAPTIONS 


Fig.l —The evolution of the density contrast (<5) with 
time for the mass perturbations in the range 10^ — 10^^ M©. 
The solid line is the internal shell, the short dashed line is the 
middle shell and the long dashed line is the external shell, 
where the time is given in years. The input parameters of 
the models are: Mo = 4 x lO^^M©, h = 1.0, Hb ~ 0.1 and 
n = —1. 

Fig.2 —The evolution of the mass density p (in g cmT^) 
for the baryonic matter with time. The captions and input 
parameters for the curves are the same as the Fig.l. 

Fig.3 —The mass density p versus radius (in pc) of the 
clouds. The solid line is the initial time (ti), the short dashed 
line, long dashed line and very long dashed line are inter¬ 
mediate times (respectively t 2 , ts and t 4 ,) and the dashed- 
pointed line is the final time (t/). The input parameters are 
the same as the Fig.l. 

Fig.4 —The evolution of the temperature (T in K), 
pressure of the gas (P in dyn cmT^) and density of molecules 
(n_ff 2 in molecules cm~^). The captions and input parame¬ 
ters for the curves are the same as the Fig.l. 

Fig.5 —The evolution of the density of molecules (riH^) 
versus temperature (T). The captions and input parameters 
for the curves are the same as the Fig.l. 

Fig.6 —The evolution of the radii of the purely baryonic 
clouds for 10^ — lO^^M© as a function of time. The solid line 
is the internal shell, the short dashed line is the middle shell 
and the long dashed line is the external shell, where the time 
is given in years. The input parameters for these models are 


© 0000 RAS, MNRAS 000, 000-000 



8 S.R. Oliveira et al. 


the same as the Fig.l. 

Fig. 7 —The evolution of the mass density (p) with time. 
The solid line is the internal shell, the short dashed line is the 
middle shell and the long dashed line is the external shell, 
where the time is given in years. The input parameters for 
these models are: Mo = 4 x lO'^^M©, h = 0.5, fls = 0.1, 
floM = 0.5 and n = —1. 

Fig.8 —The mass density (p) versus the radii of the 
clouds. The solid line is the initial time (ti), the short dashed 
line, long dashed line and very long dashed line are inter¬ 
mediate times (respectively t 2 , ts and t^) and the dashed- 
pointed line is the final time (t/). The input parameters are 
the same as the Fig.7. 

Fig.9 —The evolution of the velocity (in km s with 
time. The captions and input parameters for the curves are 
the same as the Fig.7. 

Fig.10 —The velocity profiles to the baryonic and dark 
matter components. The captions for the curves are the 
same as the Fig.8. The input parameters for the models are: 
Mo = 4x 10 ^'^Mq, h = 0.5, fisH = 0.1, Qom = 0.5 and 
n = —1. 

Fig.ll —The evolution of the temperature (T), pressure 
(P) and density of molecules (nrfs) with time of the baryonic 
matter component. The captions and input parameters for 
the curves are the same as the Fig.7. 

Fig.l2 —The temperature of the gas (T) versus the 
molecular density {uhO- The captions and input parame¬ 
ters for the curves are the same as the Fig.7. 

Fig.13 —The evolution of the radius of the baryonic 
and dark matter components of the clouds with time. The 
captions and input parameters for the curves are the same 
as the Fig.7. 


This paper has been produced using the Royal Astronomical 
Society/Blackwell Science macros. 
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